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The cause of the weak solar cycle 24 

J. Jiang^, R.H. Cameron^, M. Schiissler^ 
jiejiangSnao.cas.cn 

ABSTRACT 

The ongoing 11-year cycle of solar activity is considerably less vigorous than 
the three cycles before. It was preceded by a very deep activity minimum with 
a low polar magnetic flux, the source of the toroidal field responsible for solar 
magnetic activity in the subsequent cycle. Simulation of the evolution of the 
solar surface field shows that the weak polar fields and thus the weakness of 
the present cycle 24 are mainly caused by a number of bigger bipolar regions 
emerging at low latitudes with a ‘wrong’ (i.e., opposite to the majority for this 
cycle) orientation of their magnetic polarities in the North-South direction, which 
impaired the growth of the polar field. These regions had a particularly strong 
effect since they emerged within ±10° latitude from the solar equator. 

Subject headings: Sun: magnetic fields. Sun: activity 


1. Introduction 

Solar activity, as manifested by sunspots and bipolar magnetic regions at the solar 
surface, by energy release in flares, eruptions, coronal mass ejections, and e nergetic particles. 


varies cyclically with a (quasi-)period of about 11 years flHathawavl l2010l) . Activity cycles 


differ significantly with regard to their length and strength. The present cycle is among the 
weaker ones, particularly if compared with the sequence of strong cycles during the second 


activity (Russell et ah 201C 

i), weak polar magnetic held 

Munoz-Jaramillo et ah 2012'). and 

low heliospheric open flux (_ 

Smith & Baloeh 2008: Owens & Forsvth 

2013h. This is consistent 

with the empirical correlation of the polar held around solar activity minimum with the 


strength of the subsequent cycle fjWang &: Sheelevll2009l) . This correlation reflects the role 
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of poloidal magnetic flux connected to the polar fields as the dominant poloidal source of the 
toroi dal flux generated by differential rotation in the solar interior fjCameron fc Schiissler 
2015|), whose emergence at the surface during the subsequent 11-year cycle leads to the 
various manifestations of solar activity. 

Generally, the observed large- scale solar surface field is well represented by Surface 


Jiang et ah 


Flux Transport (S FT) simulations f Wang et al.lll989l: ISheelevll2005l: iMackav fc Yeatesll2012 


2014bl) . which address the evolution of the magnetic field as a result of the emer¬ 


gence of bipolar magnetic regions (sunspot groups) and the subsequent transport of mag¬ 
netic flux by near-surface plasma flows (convection, differential rotation, poleward meridional 
flow). Such models showed that the polar fields are reversed and built up by preferred flux 
transport of one polarity across the equator as a result of the systematic tilt of bipolar 
sunspot groups with respect to the direction of rotation. The strength of the polar fields 
mainly depends on the tilt angle, the d istribution of ernergin g flux in heliographic latitude, 
and the speed of the meridional flow (lRa.uma.nn et al.N2004ll . SFT simulations reproduce 
the polar fields (as regularly measured directly since 1976 and inferred from geomagnetic 
and solar proxies measured since the 19th century) relatively well for the cycles before cy¬ 
cle 23. However, unless they are c ontinuously adjusted by data assimilation from obser¬ 
vations flUpton fc Hathawavl l2014bl ). such simulations so far could not reprod uce the very 


weak polar field during the declining phase a nd activity minimum of cycle 23 fjjiang et ah 


flow were co 

1^ w ^ 

nsidei 

'ed as a possible remedy 

Schriiver & Liu 

2008: Wang & Sheelev 

2009; 

Nandv et ah 

2011; 

Jiang et ah 2013ll. but there is no empirical basis for changes of the kind 

required bv these models (Gonzalez Hernandez et ah 2008: 

Hathawav & Rightmire 

2010; 


Cameron fc Schiisslerl 120101) . 


While systematic variations of the tilt angles of sunspot groups ( Cameron et al.l 20101) 


were found to be an unlikely cause of the weak polar fields of cycle 23 (jSchriiver fc Liu 


2008 


Jiang et al.ll2013h . there is evidence for a considerable random element in the evolution of the 
polar fields: the observed scatter of the tilt angles causes a variation of 30-40% (standard 
deviatioi i from the meaii for many realizations) i n the resulting polar field around activity 
minima (lRa.uma.nn et al.ll2004l: Ijiang et al.ll2014aj) . Even single big bipolar regions emerging 
very near or across th e solar equator can significantly affect the amplitude of the polar field 
(jCameron et al.ll2013j) . This suggests that an answer to the question why the polar fields 
in the declining phase of cycle 23 were so weak (and, as a consequence, why cycle 24 is 
rather feeble) requires a study that includes detailed information about the individual tilt 
angles and magnetic polarities of th e bipolar regions that emerged during cycle 23. Such 
data became available only recently (iLi &: Ulrichll2012l: IStenfio fc Kosovichevll2012l) . so that, 
in contrast to previous studies, we could run SFT simulations including the actual tilt angles 
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for cycle 23. 


2. Surface flux transport simulations 


For our si i nulati ons we used the SFT code described in iRaumauu et ahl (120041 ) and 
Cameron et al.l (120 lOl) . which solves the magnetohydro dynamic indnction equation describ¬ 
ing the evolution of the radial component of the large-scale magnetic held at the solar snrface, 
B{X,(j),t) (where A and 0 are heliographic latitude and longitnde, respectively), as a resnlt 
of passive magnetic hnx transport by snrface hows, viz. 


dB 

'W 


-n{\) 


+r] 


dB 

d(j) 

1 


d 


Rq cos a dX 


[v{X)B cos A] 


d 


RX cos XdX \ dX 


L-"'© 


dB 


+ 


d'^B 


Rq cos^ a 


+ S{X, 


( 1 ) 


The latitndinal diherential rotation of Snn, hl(A), is taken from Snodgrass ( 198311 . For the 


poleward meridional how, u(A), we nsed the prohle bv ivan Ballegooiien et al.l (1199811. which is 


consi stent with the empirical prohle derived by magnetic element tracking (IHathawav fc Rightmire 
201111 . The “tnrbnlent” magnetic dihusivity of ?7 = 250 km^s“^ describes the random walk 
of the magnetic hnx elements in the large-scale convective how pattern of snpergranulation 
( Leighton 1964 1. 


The sonrce term, S{X,(j),t), represents the emergence of bipolar magnetic regions. As 
input database for these we used the sunspot observations obtained by the Solar Optical 
Observing Network (SOON) of the US Air Force, which yields 3046 bipolar regions between 
Jnne 1996 and December 2010 (2921 regions during cycle 23 until the end of 2008). The 
corresponding magnetic hnx is determined by a single parameter, i?max, which is calibrated 
by the total nnsigned snrface hux obtained from SoHO/MDI synoptic maps after rebinning 


to th e spatial resolution of the simnlation (1° in both latitnde and longitnde) (jScherrer et ah 


19951 ). The match of the observation and simnlation dnring the activity maxi r nnm o f cycle 23 
yields i?max = 592G. We applied the procedure described in lBaumann et al.l (1200411 to dehne 
the properties of the corresp onding bipolar r egion s. Further details of the source generation 
and calibration are given iu ICamerou et ahl (1201011 . 


For the tilt angles of the emerging bipolar regions we used the data set of iLi fc Ulrich 


(l2012l) . which is based upon daily magnetograms from the Monnt Wilson Observatory and 
from the SoHO/MDI. These data provide snperior information about tilt angles since they 
contain magnetic polarity information and also inclnde plage regions, both of which are not 
available when considering only white-light images of sunspot gronps (jWang et al.ll2015[l . As 
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in our previous studies, the tilt angles were mul tiplied by a factor 0.7 in order to account for 
the near-surface inflows towards active regions flGizon et alJl2010l ). 


Since single big active regions can have a signihcant impact on the polar helds fjCameron et al 


2 OI 3 I ). the input data were double-checked, thereby eliminating 22 duplicate (recurrent) ac¬ 
tive regions, which would otherwise have entered the simulation more than once - a list of 
these regions is provided in Table [H 

As a comparison case, we also carried a SFT simulation that assigns to each emerging 
bipolar region a tilt angle, a, according to a function of latitude (commonly referred to as 
Joy’s law) of the form a = T y^jXf, where the c onstant T = =1.4 and 

summation is over all bipolar regions of cycle 23 fjCameron et ahll2010l) . 


Our SFT simulations start at the end of June 1996, around the minimum of cycle 
22, using as initial condition the (interpolated and polar-held corrected) synopt ic map of 


the s urface magnetic held for Carrington rotation 1911 obtained by SoHO/MDI (jSun et ah 


201ll ). The evolution of the large-scale held resulting from hux emergence in bipolar magnetic 
regions and hux transport by surface hows according to Eq. (1) is then followed until the 
end of 2010. 


Fig. [T] shows observed and simulated time-latitude diagrams of the longitudinally av¬ 
eraged radial surface held. Although there are diherences in detail, the simulation on the 
basis of the actual individual tilt angles (middle panel) reproduce the main features of the 
evolution of the surface held rathe r well. Devia t ions i n detail can be traced back to a small 
number of errors in the tilt data of lLi fc Ulrichl f)20121 1. We have conhrmed that these errors 
do not signihcantly ahect the polar helds and dipole moment around the activity minimum, 
for which onl y the cross-equator transport of magnetic hux from low-latitude active regions 
is important f| Jiang et al.ll2ni4a|) . To avoid potential bias, we therefore chose not to correct 
these errors, which would have improved the agreement of the time-latitude diagrams. 


The main diherence between the simulation with actual tilt angles and that with Joy’s 
law becomes apparent when considering the axial dipole moment of the Sun, 


Q 

D{t) = - / (i?(A, f)) cos A sin A (JA, 

2 J-k/2 


( 2 ) 


where {B{X, t)) is the longitudinally averaged radial surface held. The axial dipole moment is 
a better suited quantity than the polar held for the comparison of observation and simulation. 
It is uniquely dehned and limits the ehect of asymmetries in the magnetic hux distribution 
on the two hemispheres (such as non-synchronous reversals of the polar helds). Around 
activity minima, the axial dipole moment is dominated by the polar helds and represents 
the relevant quantity for the poloidal magnetic hux connected to the polar helds. Likewise, 
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it also represents the open heliospheric magnetic flnx during these periods. 

Fig. [2] shows the evolution of the axial dipole moment as determined from observations 
(SoHO/MDI magnetic maps) and our SFT simulations. The simulation based on the actual 
tilt angles follows the observed evolution of the axial dipole moment well within the shaded 
error range (see below). On the other hand, employing Joy’s law for the tilt angles leads to 
a much too strong axial dipole moment in the declining and minimum phases of cycle 23. 
This behavior is typical for previous attempts to reproduce the observed weak axial dipole 
moment and polar field by SFT simulations, unless the models were arbitrarily modified 
(e.g., in terms of meridional flow variations). 


We estimate the error of the axial dipole moment in our simulations as follows. The main 
source of uncertainty are the errors in the determination of the tilt angles. We have therefore 
used a different data source, the Kitt Peak Synoptic Magnetograms, to independent l y rede - 
termine the tilt angles of about 30% of the bipolar regions considered bv iLi fc Ulrichl (1201211 . 
who used SoHO/MDI magnetograms. The regions were chosen from the SOON list, their 
area determined by applying a threshold on measured magnetic flux density, and their tilt 
then calculated using the line connecting the centers of gr a vity o f the positive and negative 
polarity parts. Comparison with the data of iLi fc Ulrichl fl2012l) yields an overall value of 
a = 11.5° for the rms scatter. However, the error is dependent on the area of the region; 
e.g., for the 50 biggest ARs, the scatter is reduced to about 5°. The dependence on area, A 
(in microhemispheres), can be roughly approximated by the function 


a{A) 


18° - 0.02A forA < 650. 
5° otherwise. 


(3) 


We carried out 20 SFT simulations with random scatter of the AR tilts according to this 
function. The rms variation of the resulting axial dipole strength is then taken as error 
estimate. 


3. Origin of the low dipole moment 


Figure 13] illustrates the reason for the discrepancy between the simulation with actual tilt 
angles and the one based upon Joy’s law. Shown are cumulative contributions of magnetic 
regions to the change of the dipole moment (since 1996) as a function of emergence latitude 
(modulus). The individual contribution of a magnetic region for la rge times is propor tional 
to A = sin(Q!) exp(—A^/110), where A is the area of the region ( Jiang et al.ll2014al) . The 
normalized cumulated contributions using the actual tilt angles (dashed curve) closely follow 
those when assuming Joy’s law (solid curve) for emergence latitudes poleward of ±10° from 
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the equator. Nearer to the equator, the curves diverge: the actual tilt angles of magnetic 
regions emerging in this range statistically deviate strongly from Joy’s law in such a way 
that the resulting total change of the dipole moment since 1996 is reduced by about 40%. 
The is mainly due to comparatively few regions in the (positive and negative) tails of the 
distribution of A, i.e., a result of random variations together with low-number statistics. 


An example for a big magnetic region with ‘wrong’ orientation (i.e., opposite to the 
majority of the magnetic regions during this cycle) of the magnetic polarities in the North- 
South direction is active region AR10696, which appeared in November 2004. Magnetic maps 
for this region are shown in Fig. HI Since it emerged near the equator, a significant amount 
of negative-polarity magnetic flux could cross the equator and thus signihca ntly reduce the 


build up of the axial dipole moment in the descending phase of the cycle fICameron et ah 


2013|). Since the magnetic flux at the poles around activity minima is typically of the order 


of the flux of only one big bipolar magnetic region, a few bipolar regions emerging near the 
equator can have a significant impact. 


4. Comparison with cycles 21 and 22 

While Joy’s law was a reasonable assumption in SFT simulations before cycle 23, we 
have seen that it is inadequate to describe the evolution of the axial dipole moment in cycle 
23. It is therefore of interest to consider a similar analysis as the one resulting in Fig. |3]for 
earlier cycles. Since magnetic information is necessary to properly determine the tilt angle 
(particularly for the big regions with ‘wrong’ tilt), we have to restrict ourselves to cycles 
21 and 22, for which Kitt Peak National Observatory (KPNO) synoptic magnetograms are 
available. Since these are based on observations along the central meridian only, they contain 
not all of the emerging bipolar regions detectable by full-disk magnetograms such as those 
from SoHO/MDI, but still provide a statistically meaningful sample. We used the synoptic 
magnetograms to determine the individual tilt angles of bipolar regions emerging in cycles 
21-23. The (normalized) cumulative distributions to the change of the dipole moment during 
these cycles are given in Figure [5l For comparison, the result for cycle 23 on the basis of tilts 
from the MDI magnetograms (Fig. [3]) is also reproduced. The two results for cycle 23 (from 
MDI and KPNO data, respectively) are consistent with each other (taking into account that 
the KPNO analysis covers only 30% of the MDI regions): both show the big deviation from 
Joy’s law for latitudes below 10°. On the other hand, the result for cycle 21 and 22 (combined 
in order to improve statistics) is fully consistent with tilt angles according to Joy’s law. This 
indicates these cycles were less strongly affected than cycle 23 by randomness concerning the 
tilt angles of low-latitude active regions in the descending phase of the cycle. 
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5. Concluding remarks 


Our SFT simulations have shown that the weak axial dipole moment (weak polar helds) 
around the activity minimum of cycle 23 are a result of a number of low-latitude bigger active 
regions whose ‘wrong’ tilt angles did not follow Joy’s law. Since the poloidal magnetic flux 
related to the axial dipole moment (directly related to the p olar helds) is the dominant so urce 
of the toroidal flux that emerges in the subsequent cycle fjCameron fc Schiisslerl l2015l) . its 
low amplitude during the declining phase of cycle 23 caused the weakness of the current cycle 
24. The dependence on the properties of big bipolar regions emerging in the vicinity of the 
solar equator, which are affected by randomness, strongly limits the scope of predicting solar 
cycle strength: until the amplitude of the polar fields is established around solar minimum, 
no sensible prediction of the strength of the next cycle can be made. Sound predictions 
projecting even further into the future are virtually impossible. 


The authors are grateful to Jing Li for providing the tilt angle data in electronic form. 
SoHO is a project of international cooperation between ESA and NASA. The SOON data 
were downloaded from http://solarscience.msfc.nasa.gov/greenwch.shtml, the SoHO/MDI 
synoptic maps from http://solarscience.msfc.nasa.gov/greenwch.shtml, and the KPNO syn¬ 
optic maps from http://dighb.nso.edu/ftp.html This work utilizes SOLIS data obtained by 
the NSO Integrated Synoptic Program (NISP), managed by the National Solar Observatory, 
which is operated by the Association of Universities for Research in Astronomy (AURA), 
Inc. under a cooperative agreement with the National Science Foundation. 

J. Jiang acknowledges support by National Science Foundation of China (grants 11173033, 
11221063, 2011CB811401) and support by the Strategic Priority Research Program of the 
Chinese Academy of Sciences (grant XDB09040200). 
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Table 1: List of the recurrent ARs. At Time (yyyymmdd), the ARs reached their maximum 
area (in microhemispheres). The column AR gives the active regions that are retained in 
the SFT simulations._ 


No. 

Recurrent AR 

AR 

Time 

Latitude 

Area 

Tilt 

1 

8375 

8398 

19981109 

21.0 

964 

48.13 

2 

8771 

8739 

19991125 

-15.0 

1097 

3.63 

3 

8759 

8731 

19991109 

10.0 

1230 

1.54 

4 

9085 

9046 

20000720 

15.0 

405 

29.88 

5 

9131 

9090 

20000825 

15.0 

312 

38.84 

6 

9636 

9601 

20010926 

14.0 

538 

36.56 

7 

9670 

9628 

20011019 

-18.0 

711 

16.66 

8 

9672 

9704 

20011026 

-18.0 

791 

44.88 

9 

9684 

9715 

20011104 

5.0 

737 

16.04 

10 

9751 

9715 

20011226 

4.0 

671 

15.86 

11 

9973 

9934 

20020529 

-16.0 

1283 

8.43 

12 

10036 

10069 

20020721 

-6.0 

1429 

17.57 

13 

10105 

10069 

20020910 

-8.0 

2028 

28.38 

14 

10386 

10365 

20030617 

-7.0 

418 

54.31 

15 

10464 

10484 

20030926 

4.0 

817 

17.88 

16 

10507 

10488 

20031119 

10.0 

1190 

13.23 

17 

10508 

10486 

20031119 

-17.0 

937 

25.08 

18 

10520 

10484 

20031216 

1.0 

232 

56.46 

19 

10661 

10652 

20040818 

10.0 

724 

24.04 

20 

10708 

10696 

20041203 

8.0 

219 

55.61 

21 

10743 

10735 

20050319 

-8.0 

525 

40.70 

22 

10750 

10735 

20050408 

-7.0 

205 

51.75 

23 

10775 

10792 

20050612 

10.0 

485 

31.27 

24 

10803 

10792 

20050825 

12.0 

272 

40.14 

25 

10875 

10865 

20060425 

-10.0 

644 

32.90 

26 

10904 

10908 

20060812 

-13.0 

1003 

9.94 

27 

10930 

10923 

20061212 

-6.0 

910 

29.36 

28 

10935 

10923 

20070104 

-6.0 

365 

41.59 
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Fig. 1.— Time-latitude diagrams of the longitudinally averaged radial magnetic field at the 
solar surface. Top panel; observation (SoHO/MDI synoptic magnetic field data); middle 
panel; simulation using the actual tilt angles of bipolar magnetic regions; bottom panel; 
simulation using tilt angles according to a htted latitude dependence. 
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Fig. 2.— Time evolution of the solar axial dipole moment. The curves correspond, respec¬ 
tively, to observed SoHO/MDI magnetic maps (black), a simulation using the actual tilt 
angles of bipolar magnetic regions (red), and a simulation using tilt angles according to a 
htted latitude dependence (blue). 
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Fig. 3.— Cumulative contributions (in latitudinal distance from the equator) of magnetic 
regions to the change of the axial dipole moment since 1996. Results are given assuming 
tilt angles according to Joy’s law (oc solid curve) and for the actual individual tilt 

angles (dashed curve). 
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Fig. 4.— Example of a bipolar magnetic region that significantly weakened the axial dipole 
moment in the declining phase of cycle 23. Shown are SoHO/MDI magnetic maps of the 
active region AR10696 taken Nov 5, 2004 (left panel) and Dec 2, 2004 (right panel, after one 
solar rotation, then denominated AR10708), respectively. Positive magnetic flux is indicated 
in white, negative flux in black. Owing to its near-equator emergence, high tilt, and abnormal 
polarity orientation in the North/South direction, the region provides a signihcant amount 
of negative flux that is transported over the equator (by supergranular random walk) to the 
southern hemisphere. Through poleward advection (mainly by meridional flow) this flux 
eventually weakened the buildup of positive flux around the south pole of the Sun, thus 
lowering the axial dipole moment. 
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Fig. 5.— Cumulative contributions (in latitudinal distance from the equator) of magnetic 
regions to the change of the axial dipole moment during cycles 21, 22, and 23. Full lines 
give the result assuming tilt angles according to Joy’s law (almost coincident with each 
other), dashed lines for the case using the actual individual tilt angles. The black curves are 
reproduced from Fig. [3] (on the basis of MDI full-disk magnetograms), blue lines indicates 
results from KPNO synoptic magnetograms for cycle 23, and red lines the corresponding 
results for cycles 21 and 22 combined. 





